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Abstract — This paper deals with the accomplishment of total 
area coverage of an arbitrary region using sensors with a finite 
sensing radius of r^. For a given region, we aim to obtain a 
deterministic placement of sensors which, apart from ensuring 
that the entire region comes under the purview of at least 
a single sensor, minimises the number of sensors utilised. We 
begin by considering regions devoid of obstacles and thus having 
every location amenable for placement. Herein, we formalise the 
popular notion that sensors at the centres of the hexagons of 
a hexagonal tessellation provide the most optimal placement. 
We then move on to regions which may comprise obstacles of 
arbitrary size at arbitrary locations. We recognise two distinct 
classes of obstacles, namely transparent and opaque obstacles, 
which are distinguished by their ability (or the lack of it) 
to permit sensing radiation through them. In the real world, 
transparent obstacles model lakes, ponds and swamps, while the 
opaque ones stand for, inter alia, hills, trees and walls. We propose 
a polynomial-time algorithm for achieving optimal placement in 
the aforesaid scenarios and we prove its convergence. 

Index Terms — Coverage, Transparent and opaque obstacles. 
Art gallery problem 

I. Introduction 

A wireless sensor network (WSN) consists of a large 
number of tiny, low-powered devices, called sensors, which 
communicate with each other (possibly) in a multi-hop fash- 
ion. These sensors suffer from severe constraints, such as low 
reliability and limited battery power (or energy), processing 
power, storage memory, sensing range and communication 
range capabilities, to name a few, with energy being the most 
critical one. Tremendous advances in inexpensive sensor tech- 
nology and wireless communications have rendered the design 
and development of large-scale wireless sensor networks cost- 
effective and viable enough to attract the attention of a wide 
range of civilian, natural, and military applications. 

One of the fundamental issues that arises in sensor networks, 
in addition to location calculation, tracking and deployment, 
is coverage. Coverage can be construed to be a measure of 
quality of service of a sensor network, in the sense that a 
high degree of coverage of a region signifies that the region is 
monitored well. Let us consider a specific application where 
coverage of a region becomes a critical necessity. Millions 
of acres of land are lost every year around the world due to 
forest fires. An early detection of these fires plays a crucial 
role in preventing them from causing irremediable calamity. 



It is possible to employ specialised sensor nodes in critical 
areas to monitor them and thereby preclude the occurrence 
of such disasters. The benefits of the placement of a large 
number of nodes to cover such areas clearly override the cost 
factor It is fairly evident that coverage is a QOS parameter 
in this example, as the odds of a fire being detected early is 
proportional to the extent of coverage employed. 

The problem of area coverage is thus elementary in the 
design of wireless sensor networks. In this paper, we aspire to 
design algorithms that could achieve total area coverage in an 
arbitrary two-dimensional region. A two-dimensional region is 
said to be 1 -covered if every point (excluding the obstacles) 
inside the region is within a distance of Vg from its nearest 
sensor, where is the sensing radius of the sensor Wireless 
sensor nodes (or motes) have to be placed in the region to 
obtain complete coverage. The sensor design and the sensing 
modality is such that the coverage region of a sensor is a 
disk of radius centred at the location where the sensor is 
placed. The region to be covered has transparent obstacles 
(e.g., lakes, ponds and swamps) where the sensors cannot be 
placed but through which the signal being sensed can pass, 
as well as opaque obstacles (e.g., trees and other structures), 
which neither allow the sensors to be placed nor allow the 
signal being sensed to penetrate through them. The problem 
we deal with is to cover the area not occupied by obstacles 
with as few as sensors as possible, ensuring that a specific 
coverage objective (e.g., fc-coverage) is met. 

II. Related Work 

There are three main classes of coverage problems tackled 
in the literature on coverage. These are barrier coverage, area 
coverage and point coverage. Barrier coverage involves placing 
sensors along vulnerable boundaries of a region for achieving 
intrusion detection. Area coverage refers to the placement 
of sensors to monitor the entire area, while point coverage 
focusses on the monitoring of certain critical points identified 
within a region. We now cast a quick glance over some of the 
specific coverage problems dealt with in the literature. 

Gautam Das et al. |10| provide a O(fclogfc) algorithm to 
evaluate the placement of an optimal number of sensors to 
cover a given region without obstacles. An initial random 
placement of the k sensors is followed up with the Voronoi 



partition of the area. In every iteration of the algorithm, the 
position of the sensor is moved towards the circumcircle of 
the Voronoi polygon. It is proved that the new position of 
the sensor lies within the same polygon and also that the 
radius of the circumcircle reduces with every iteration. The 
procedure terminates when the radius of the circle does not 
shrink anymore. 

The above algorithm can be adapted to solve the problem 
of optimal sensor placement in a region without obstacles. We 
can start with an arbitrary number of sensors and follow the 
procedure. An appropriate value of the number would be the 
ratio of the area of the region to the sensing area of a sensor 
If the algorithm terminates with the maximum radius of all 
the circumcircles being greater than the typical sensing radius 
of a sensor, the iteration needs to restart with more number of 
sensors. 

Senjuti Basu Roy f^'l considers in detail the problem of 
finding the best coverage path in the presence of transparent 
and opaque obstacles. The best coverage path (also called the 
Maximal Support Path) is defined as the path between a source 
and a destination point which has the minimum cover value, 
where cover value of a path is obtained by maximising over the 
distances of the points on the path from their corresponding 
nearest sensors. For the case of opaque obstacles, the crux 
of the algorithm presented is to compute the dual of the 
Constrained Voronoi Diagram (CVD) and to run a Bellman- 
Ford algorithm along the edges of the dual to find the shortest 
path. In the case of transparent obstacles, a visibility graph is 
used in place of the CVD. 

Murray et al. tSJ extend the VDH-based solution [lO] to 
non-convex regions. The 1 -centre problem, which is funda- 
mental to solving the above problem, has been solved using 
efficient algorithms. The p-centre problem has the closest 
resemblance to the problem of total area coverage of a region 
A. The p-centre problem is one of optimally placing the p 
facilities (or sensors, equivalently) so that the distance of any 
point in A to its nearest sensor is minimised. It has both 
discrete and continuous versions. In the discrete version, the 
task is to select a subset of p points from a large set such that 
placement of facilities on the points of the subset minimises 
the maximum distance of any other point in the set to its 
closest facility. The continuous version relaxes the constraint 
of having to place the facilities at the demand points. Instead, 
an underlying graph or map is defined and the facilities can be 
placed anywhere along the edges of the graph or the regions 
in the map so as to minimise the maximum distance of any 
demand point to its nearest facility. 

The continuous p-centre problem may be formally stated as 
follows: 



min max [ mmdj{x,y) 

(xj,yj)eA, j=l,2...,p \{x,y)eA\ ] j 



where 



A : the analysis region 

p : no. of facilities located 

{x, y) : reference point in A 

In terms of the p-centre problem, the problem of covering a 
plane region without obstacles reduces to 



mm p 

C<r^ 



where 



C :— min 



max mm dj{x,y) 



where is the sensing radius of a sensor 

The continuous p-centre problem is proven to be NP-hard 
Yet the minimisation problem of the number of sensors 
with respect to the radius constraint has a surprisingly simple 
solution, as we will prove in the ensuing section. 

III. Coverage of a Plane Region sans Obstacles 

The problem of finding the minimum number of circles to 
cover an arbitrary two-dimensional region is well-studied and 
it is known that a tessellation of regular hexagons (inscribed 
in the coverage circle) solves the problem. Kershner [l], in 
1939, gave the first formal proof of such a result, when the 
radius of the coverage circle reduces to zero for a fixed area 
to be covered. Particularly, he proved the following result. 

Theorem III.l. Let M denote a bounded and closed plane 
point set, N{r) be the minimum number of circles of radius 
r which can cover M and A{M) represent the area of M. 
Then, 

lim TTr^N{r) = ^^A(A/) 

7 — S-0 9 

We note that the minimum number of disks of radius r 
required to cover M is "^^^P . However a certain amount 
of overlap between neighbouring circles cannot be avoided, 
thus necessitating more disks than this minimum. Kershner's 
theorem actually says that the number of disks required is 
greater than this minimum by a factor '^^^ {— ^ ), which 

is the ratio of the area of a circle of radius r to its inscribed 
hexagon. We extend this result to establish the optimality of 
hexagonal coverage of a rectangle using circles of a constant 
radius, as the area of the rectangle grows to oo. 

A. Optimality of Hexagonal Tessellation for Infinite Plane 
Areas 

Theorem III.2. Let I and w be the length and width of a 
rectangular region A of area A (A) be the number of 

circles of radius rg that cover A and Aq = ^^^s- Then, 

i->-oo -p- 

Proof: The result follows from Lemmas 5 and 6 in |[T1. 
From Lemma 5, 



dj{x,y) : distance of point (x,y) to the j facility 
{xj,yj) : coordinates of facility j 
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From Lemma 6, we have, for perimeter p of A, 

irrlN}^ [A) < ^ {A + 2pr, + IGr^) 
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which completes the proof. ■ 
Similar bounds can be established for coverage of A with 
hexagons. 

Theorem III.3. Let I and w be the length and width of a 
rectangular region A of area A, Nq (A) be the number of 
hexagons of circumradius that tessellate A and Aq = 
2^r2. Then, 




Fig. 1 . Peripheral Hexagons displaying the finite edge effects which establish 

Nff (A) 

the lower bound for — ^^ 



Proof: The finite edge lengths dictate that not all 
hexagons employed of the tessellation occupy an area of Aq. 



NSiA) 



> 1 



Np;(A) 
lim inf ".^ ' > 1 



Let fci and ko be constants such that -r-^ and refer 
to the number of hexagons along the length and width of 



the rectangle respectively, fci and k2 are dependent upon the 
position and orientation of the tessellation used. 

/ in 

iVg(A)v4o - 2^0(7 + ) < A 



'kiTs k2rs 

,. Nl^jA) 
lim sup — 



< 1 



Corollary III.4. For a two-dimensional infinite plane region, 
tessellation with hexagons is optimal with regard to the 
number of hexagons used, i.e.. 



If;— ^00 w 
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IV. Coverage in the Presence of 
Transparent Obstacles 

In this section, we concentrate on extending Kershner's 
work to regions with transparent obstacles at arbitrary loca- 
tions and of arbitrary shapes and sizes. 

A. Terminology 

Kershner's result establishes the fact that for an infinite 
plane area sans obstacles, tessellation with hexagons provides 
an asymptotic minimum with regard to the number of sensors 
used. Such a universal hexagonal tessellation (a hexagonal 
tessellation spanning the entire 2-D plane) is defined by 
two parameters: the centre of one of its hexagons and the 
orientation of the hexagons (i.e., say the angle made by any 
one of the diameters of the hexagon with the positive x-axis.) 
In this section, let U denote one such universal hexagonal 
tessellation of an infinite plane area. An area is said to be 
occupied by a transparent obstacle if it is unavailable for 
placement of sensors, but is not impervious to the signal 
to be sensed. Let be a finite and bounded plane area 
interspersed with transparent obstacles. Obstacle-free regions 
within A are referred to as 'lands' and are accessible, while 
those regions that are occupied by transparent obstacles are 
termed inaccessible. An island is defined in the usual sense, 
i.e., an area of land surrounded by a transparent obstacle all 
around it. A(^) stands for the accessible area present within 
the region A. 

B. Bounds on the Number of Sensors 

A is initially tessellated with hexagons from U, and since A 
has transparent obstacles, some of the hexagons' centres may 
be found to lie in inaccessible regions. We call this initial 
tessellation the primary tessellation. Some of the hexagons 
may lie entirely within the transparent obstacle and hence 
the corresponding sensor is not required to cover any land. 
Hexagons whose centres lie in a transparent obstacle, but 
which have some land requiring to be covered are termed 
anomalous hexagons. Let A denote the sum of areas of all 
normal hexagons whose centres are accessible or 'intact', and 
Ao denote the sum of areas of all anomalous hexagons. The 



following theorem provides us with bounds on the number of 
sensors. 

Theorem IV.l. For an area A with transparent obstacles, 
which has been tessellated with lA, the number of sensors, 
N, required for achieving total area coverage is bounded as 
follows: 

A + 5A„ 

N < —-. (2) 



N > 



A + Aq 
Ao 



(3) 



where 

A : sum of areas of all normal hexagons 
Ao : sum of areas of all anomalous hexagons 

An intuition about the extreme situation can be gained from 
Figure |2] In the anomalous hexagon shown, we identify five 
five possible land locations and assume that the rest of the 
hexagon is occupied by the transparent obstacle. The circles 
shown are the coverage disks of the sensors placed at these 
locations. It could be seen that none of the five sensors aid 
in covering more than one land location, thereby making their 
usage indispensable. We need the following lemmas in order 




• Potential Land Locations 

Fig. 2. Worst Case Positions of Islands witliin a Hexagon 

to prove Theorem II V. 1 1 

Lemma IV.2. No six points can be located in the closure of 
a regular hexagon of circumradius such that every pair of 
the points is at a distance > from one another. 

Proof: We prove the lemma by contradiction. Suppose 
there exist six such points xi...xq. They clearly belong 
to different sextants within the hexagon, as the maximum 
distance between two points in a sextant is r^. Identify the 
sextants to which they belong. Drop a perpendicular from 
each of these points to the edge of the hexagon defining its 
corresponding sextant. 

Claim: If we move the points xi . . .xq along the perpen- 
diculars to meet the edges at say, x'l . . . Xg, distance between 
no pair of points decreases. 

Proof of Claim: 
We identify three cases for the angles with which the perpen- 
diculars meet: 




Fig. 3. Tliree cases for the points lying inside different sextants 



• The perpendiculars are either parallel or meet at 180°, in 
which case the claim is obvious. 

• The perpendiculars meet at 120°. If a, b are the initial 
distances of the points from the intersection of the perpen- 
diculars and e is the distance moved by the points along 
their respective perpendiculars (e > 0), from cosine law, 

+ + ab<{a + ef + {b + ef + [a + t){b + e) 

• The perpendiculars meet at 60°. If the points, say xi and 
X2 are at distances a and b respectively from the point of 
intersection, the distance between them is + 6^ — ab. 
So, 

+ (a + 6)e > ^ 

+ b^ -ab<{a + ef + (5 + ef - (a + e)(6 + e) 

Thus the claim is proved in each of the three cases and 
hence the six points could be safely assumed to be on the 
sides of the hexagon. Clearly, each of these points must lie 
on a different edge. Now join the six points to form a convex 
polygon. The polygon is convex since all interior angles are 
< 180°. 

Claim: If a convex polygon lies in the closure of another 
convex polygon, the inner one has a smaller perimeter than 
the outer one. 

Proof of Claim: Consider the six triangles, each of which 
shares a unique edge with the interior polygon, and has the 
other two edges along the perimeter of the outer polygon. 
Triangle law can be invoked in each of those triangles and the 
resultant inequalities can be added up to prove the claim. 

The perimeter of the assumed convex polygon is > Gr^, 
which is a contradiction to the previous claim. The existence 
of such a polygon is thus invalidated. ■ 

Lemma IV.3. Five points can be located in the closure of a 
regular hexagon of circumradius Vs such that no two of them 
are at a distance < from each other 

Proof: The proof follows from basic geometry and is 
omitted for the sake of brevity. ■ 

Proof of Theorem \IV.1\ The bounds can now be estab- 
lished using the previous Lemmas. We just need to observe 
that the hexagons under Ao necessitate at most five sensors 
each. The land areas that make up A require just one sensor per 
hexagon (Figure|2]i. Similarly, a lower bound on the number of 
sensors for covering any region arises in a land pattern where 
the anomalous hexagons demand no more than a single sensor. 



C. 1-Coverage of Anomalous Hexagons 

We propose an algorithm for coverage in the presence 
of transparent obstacles that involves recursively sliding the 
universal hexagonal tessellation U over the area that remains 
to be covered iteratively until no region is left uncovered. We 
recall the earlier notation of A for sum of areas of hexagons 
with centre intact, i.e., normal hexagons and Ao for sum of 
areas of hexagons with centre inaccessible but some land left 
to be covered, i.e., anomalous hexagons. 

We define a cluster as a group of anomalous hexagons 
contiguous to one another, where two hexagons are said to 
be contiguous if they share an edge. Let Pi, P2, ■ ■ ■ Pk be 
the clusters formed among the ^g- anomalous hexagons. It 
could turn out that some (or even all) of the clusters may 
consist of a single hexagon. The goal now is to reduce 
the uncovered area in each of those clusters. The algorithm 
works separately on each cluster to reduce its uncovered 
area. Consider the following seemingly obvious fact about a 
hexagonal tessellation: When the centres of all the hexagons 
comprising a universal hexagonal tessellation are shifted by a 
specific distance and angle to points within their corresponding 
hexagons, the resulting points can be interpreted to be the 
centres of a new universal hexagonal tessellation. The common 
feature in all anomalous hexagons is the inaccessibility of their 
centres. This property enables us to find a suitable shift to their 
centres such that the maximum area is covered. 

Lemma IV.4. Consider a hexagonal tessellation o/R^. If the 
centres of each of the hexagons is shifted by a distance and 
angle to another point within the same hexagon, these shifted 
centres form a new hexagonal tessellation of K.^. 

D. Analysis of the Algorithm 

The convergence of the algorithm can be established by 
noting that in every iteration, the superposed hexagonal tessel- 
lation will reduce some area that is yet to be covered, though 
perhaps not all. Let us denote by L the set of all possible shifts 
of the centre of a hexagon to another point within the same 
hexagon. We note that L is uncountably infinite and a shift is 
defined by the distance and angle with respect to a reference. 
Let I S L represent one such shift. Let Ui : K+ — > M+ be 
a function that operates on a region with a specific area. If 
A represents the uncovered region in a cluster of anomalous 
hexagons, Ui{l^{A)) denotes the area remaining to be covered 
in A after the sensors at the centres of the anomalous hexagons 
in the cluster are shifted by I. Three properties of Ui are 
identified: 

1) ¥ov A<Z A and ^{A) ^ 0, Ui{^{A)) < A{A) with 
strict inequality holding for at least one I. 

2) UiiA{A)) = A(i), V / ^ A{A) = 0. 

3) If Ac A, A' d A and An A' = 0, Ui{A{AuA')) = 
UiiAiA)) + Ui{A{A')) 

The first property holds because as long as there is some 
uncovered region in any of the hexagons, the centres could 
be shifted to a point within them thus ensuring a strict 
reduction in uncovered area. The second statement follows 



as a corollary, because if no shifting of centres yields any 
reduction in uncovered area, it means that the entire region 
is occupied by a transparent obstacle. The third and obvious 
fact justifies dealing with clusters of anomalous hexagons. 
The area reduction achieved by a particular shift of centres 
of geographically separated clusters of hexagons is equal to 
the sum of reductions achieved by shifting each of them 
separately. 

The algorithm works as follows. Every cluster of hexagons 
is considered separately. For every cluster, we choose a par- 
ticular shift of the centres of hexagons from the set L that 
leaves the minimum area uncovered. We understand that it is 
an infinite minimisation requiring approximation in order to 
be executed in real time. One iteration involves computing 
the optimal shifts for all clusters in the region and locating 
the sensors. Some of the anomalous hexagons may have been 
covered completely during the iteration. We neglect them and 
we identify the clusters among the remaining ones. We follow 
the same procedure until the entire region is covered. We 
now analytically prove that this algorithm stops after a finite 
number of iterations. 




Anomalous Hexagons 
Normal Hexagons 
Transparent Obstacles 



Areas uncovered after 1 st iteration 
Node locations of U in 1 st iteration 



(a) Initial location of obstacles (b) After an Iteration of the Algo- 
rithm 

Fig. 4. Depiction of the Algorithm 



We define A{R^) to be the area of the region remaining to 
be covered after the j*'* iteration and A(i?^) as the area of the 
region that remains to be covered in the i*'' cluster after the 



J 



th 



iteration. Let TV^^' be the number of clusters remaining 



after the j iteration. We have, 



A{W) = A{W-^)- 



E 

1=1 



(^A(i?r')-minC/,(A(i?r') 



where 

i=l 

A(i?^H E inmC/,(A(i?r')) 

i=l 

^0-1) 

<min E C^^(A(i?r')) 

i=l 

The penultimate equation tells us that the minimum uncovered 
area obtained by shifting the centres of anomalous hexagons 
spread over the entire cluster is greater than the minimum 
obtained by considering individual clusters. The last equation 
comes from the definition of a cluster and the third property. 
Therefore, 



A(i?^') < Ui{A{R^-^)) V I 
< A(i?^-i) 



(4) 
(5) 



for at least one I, if A(i?^ ^) 7^ 0. This follows from the 
first property. 

Lemma IV.5. The algorithm stops in a finite number of 
iterations if the number of anomalous hexagons is finite. 

Proof: The algorithm achieves strict reduction in area 
remaining to be covered in every iteration. This means that at 
least one sensor is placed in the course of an iteration. But the 
maximum number of sensors required in any area is bounded 
by Equation 121 Hence the number of iterations to convergence 
is also bounded by ■ 

Lemma IV.6. The sequence {A(i?-')} converges to zero. 

Proof The sequence {A(_R^)} is a strictly decreasing 
sequence all of whose members are finite and positive. We use 
Lemma HV.SI and the algorithm stops in N iterations. Suppose 
say, for contradiction, the sequence has a positive limit p > 
and that further iterations of the algorithm do not cause a 
reduction in uncovered area. 

lim A(W) = p 

A(i?^+i) < Ui{A{R^)) V I 

< Ui{p)<p 

Ui{p) < pandUi{p)>py I 

Ui{p) = p\f l^ p = 0. 

The second equation follows from Equation |4l whereas the 
third is a result of the first property of [/;. The last statement 
follows from the second property of [/;. ■ 

E. Pseudocode and Complexity 

In practice, the optimisation can be carried out over only 
a finite number of shifts. Hence we define a lattice of points 
inside a hexagon which represent the results of permissible 



shifts of the hexagon's centre. Thus we optimise over the set 
L' C L. We begin with an initial hexagonal tessellation to 
cover the entire area inclusive of transparent obstacles. We 
then classify the hexagons as normal or anomalous depending 
on whether their centres are accessible or not. 

Let A be uncovered region in a single cluster and A be 
its area. Let AI be the number of clusters after the initial 
hexagonal placement and \L'\ be the number of lattice points 
within each hexagon. Let Ui{A) be the area remaining to 
be covered in a region with area A after placement of a 
tessellation the centres of which are obtained by applying a 
shift I E L' on the centres of the anomalous hexagons. Let 
us define X to be the set of all coordinates of centres of 
hexagons when a shift I E L' is applied to the existing centres. 
Let X be an element of X and H{x) refer to the hexagon of 
circumradius centred at x. The following equation computes 
UiiA). 



Ui{A) = A-A 



U(^(^)W}))n-^ 



X 



(6) 



In this equation, the second term on the RHS represents the 
area that gets covered because of the new tessellation. Some 
of the centres of the new tessellation could also be lying in the 
transparent obstacle and they may not contribute to coverage. 
The areas of those hexagons are excluded using the expres- 
sion l{xeA}- Also the areas of some of the new hexagons 
could partially overlap with regions that have already been 
covered. These are excluded from contributing to coverage 
again by means of the intersection with A. The following 
pseudocode summarises the minimisation procedure adopted 
at every cluster The result of the pseudocode, best_shift 
gives the particular shift of hexagons' centres that yields the 
maximum reduction in uncovered area in the cluster. 
Pseudocode : 

Ubest_shift{A) = 00 

for i = 1 to \L'\ 

hest_shift=aig v[\\iv{Uhest_shi ft {A), Ui{A)) 

end 

The above procedure is repeated at each of the M clus- 
ters present at the beginning of a particular iteration. If the 
number of anomalous hexagons in a cluster is given by iV^, 
i = 1,2... M, the complexity of the algorithm on a cluster 
i is 0(^Ni\L'\'^y The overall complexity per iteration then 
becomes 0((^^^ A^i)|L'p), where J2iLi gives the total 
number of anomalous hexagons distributed over the M clus- 
ters. The overall complexity of the procedure, given it stops in 
N iterations, is then 0((I]f=i Ni)\L'\'^N). Strictly speaking, 
the number of anomalous hexagons in every iteration reduces; 
yet the expression serves as an upper bound. If we take the 
number of permissible shifts of the centre of a hexagon to 
be a constant, the complexity becomes 0((^*^j^ Ni)N'^ = 
O(-r^N), where -4^ gives the total number of anomalous 
hexagons present at the start of the procedure. 



V. Coverage in the presence of Opaque Obstacles 

A general version of problem of coverage in the presence 
of opaque obstacles has been studied previously in literature 
under the tag of art gallery algorithms. The naive art gallery 
problem is posed as follows: 

'If a location is said to be monitored by a particular 
guard if a direct line-of-sight path exists from the guard to 
the location, what is the minimum number of guards required 
to monitor a given polygon? What are their corresponding 
placements within the polygon?' 

The original art gallery problem dates back to 1973 and to- 
day the literature abounds with research articles on algorithms 
to solve the original problem and its myriad variants. The 
opaque obstacle problem can be mapped to one such variant - 
the art gallery problem with bounded visibility (i.e., since the 
sensor has a given coverage radius). The extra constraint of 
the guards having finite visibility is added to the above stated 
problem. 

A polynomial time algorithm to this end was proposed 
in lfT2]| . The algorithm essentially involves convexification of 
the input polygon and identifying a suitable location within 
each convex polygon to place a guard. The main downside 
of the algorithm is that some of the convex polygons tend 
to be extremely small and the utility of the guards present 
in such polygons is reduced. Moreover, the union of a few 
adjacent polygons may be monitored by a single guard in some 
circumstances, though the union may not necessarily remain 
convex. Bhattacharya et al. in lfT4ll propose an incremental 
algorithm for coverage. The algorithm starts with an arbitrary 
point and computes the visibility polygon of the point. The 
next point is chosen to be on the boundary of the first point. 
The procedure is repeated until the entire region is covered. 
It is evident that this algorithm also trades off the number of 
sensors for guaranteed coverage. 

A. Problem Formulation 

We borrow two terms from the jargon of art gallery liter- 
ature. A star polygon is a polygon which can be monitored 
using a single guard, i.e., there exists at least one location 
inside the polygon such that a line-of-sight path can be 
constructed from that point to every other point in the polygon. 
Such a point is called a kernel. It may be noted that all convex 
polygons are star polygons and all points in its interior are their 
kernels. Let A denote the region that requires sensor coverage. 
We define the 'restricted star polygon (RSP)' of a point x as 
follows: 

RSP{x) = {y €A\ \\y - x\\ < r„ Xx + (1 - X)y € A 

VAe [0,1]} 

It is to be noted that the term 'polygon' is a misnomer. The 
RSPs are circles for points on an infinite plane region without 
opaque obstacles. Figure |5] depicts the RSP of a sensor in 
the presence of an opaque obstacle. In terms of the RSP, the 
problem at hand can be expressed as follows: 

min iV (7) 




• Sensor Norte 
^1 Obstacle 

Ai^^^ liiddoii by obstacles 
I I Visible aren 

Fig. 5. RSP of a sensor in tlie presence of an obstacle 



The goal as stated by the above formulation is to find the 
minimum set of points, the union of whose RSPs span the 
entire area A. 

B. Complexity 

The general art gallery problem and many of its variants 
have been established to be NP-hard. In this section, we prove 
the NP-hardness of the stated problem. 

Theorem V.l. The problem of total area coverage of a region 
with opaque obstacles is NP-hard. 

Proof The proof involves a mapping to the minimum set 
cover problem. The minimum set cover problem can be posed 
as below: 

'Given a finite set S, a collection of subsets of S, namely 

C, and an integer N, does there exist a C such that C C C, 
UC" = S and \C'\ < N?' 

Let A be the region requiring coverage. To cast the decision 
version of the opaque obstacle problem as a min-set cover 
problem, we identify the following mappings: 

A S 
RSP{x),x E A sets constituting C 

Are there N points such that A C u'jLj^RSP{xi)? ^ 
Min Set-cover problem 

The opaque obstacle problem is certainly in NP, as a 
polynomial-time algorithm could be devised to check if the 
RSPs in the solution combine to give A. Hence the decision 
version of the problem is NP-complete and the optimisation 
problem is NP-hard. ■ 

C. Bounds on the Number of Sensors 

The bounds on the number of sensors required by a region 
containing opaque obstacles do not exist in the absence of 
some tailor-made assumptions. This can be intuited by ob- 
serving that for every sensor placed inside the hexagon, the 
coverage zone of the sensor can be limited to an infinitesimal 
area around it by introducing an opaque obstacle that engulfs 
it on all directions. Three assumptions are made to circumvent 
the above stated hitch in the derivation of bounds. Let Hn^^ (x) 



denote a hexagon centred at x and of a given circumradius 
Ro, Ro << rg. The bounds on the number of sensors can be 
obtained in terms of Rq- Let K represent the obstacle. 

• if is a convex set. 

• The obstacle must be large enough so that a hexagon 
of circumradius Ro fits inside it. In other words, the 
obstacle K is deemed valid if and only if there exists at 
least one x £ K such that Hj^^ (x) C K. 

• None of the sensors employed should cater to the cov- 
erage needs of a region smaller than a hexagon of 
circumradius Rq- In other words, we do not care about 
the coverage of a region that does not contain at least 
iiflo . Consider an anomalous hexagon where the entire 
region is occupied by opaque obstacles, except for a mass 
of land A„i. The assumption is that there exists at least 
one X E Am such that Hr^{x) C Am- 

The bounds can be derived in the light of the above assump- 
tions. Let A be the region for which total area coverage is 
desired. We begin by tessellating A with an arbitrary hexago- 
nal pattern U. Complementary to the transparent obstacle case, 
we proceed to define 'anomalous hexagons'. 

Definition V.2. Ifx is the centre of one of the hexagons formed 
by the tessellation, the hexagon Hr^ (x) centred at x and of 
circumradius rg is said to be 'anomalous' if either of the 
following conditions hold. 

. xe Abut RSP{x) n Hr^ {x) ^ Br, (x) n A 

• X ^ A and 3 at least one y, y £ H^, {x) and y £ A 




Area coveredby a central Kex^n 

■ Area imcoveredby a central liex^n 

■ Opaque otstaele 



Fig. 6. Anomalous hexagons 

The first proviso is applicable whenever we have a hexagon 
whose centre is on land, but has some opaque obstacle lying 
elsewhere within it. The second condition applies to hexagons 
whose centres are on obstacles, but there is some region 
uncovered within the hexagon. The area of such regions are 
lower bounded by our assumptions. Figure |6] presents three 
different instances of anomalous hexagons, the first of which 
falls under the first alternative and the rest under the second 
alternative of the definition. 

Tlieorem V.3. For a region A with opaque obstacles, which 
has been tessellated with U, the number of sensors, N, for 
total area coverage is bounded as follows: 



- Ao 

where n is the largest number of hexagons of circumradius 
Ro that could be present (partially or completely) inside a 



hexagon of circumradius Vg. 



- Ao 

Proof: Consider a single anomalous hexagon. We need 
to compute the maximum number of sensors that may be 
required by this hexagon. We create a tessellation of hexagons 
of circumradius Ro within the hexagon. Let the total number 
of such hexagons be n. The maximum number of sensors 
necessitated by the rs -hexagon can be computed by identifying 
the corresponding distribution of obstacles and plane areas 
within the hexagon. The assumptions that we have made 
ensure that the smallest area that a sensor may be required to 
cover is equal to the area of a hexagon of circumradius Ro. 
This happens when one of the hexagons in the tessellation 
has uncovered area, but all of its neighbours consist of only 
obstacles. The maximum number of sensors are demanded 
by a region whose uncovered area appears as hexagons of 
circumradius Ro and are duly covered by obstacles all around. 

Thus in the tessellation, we segregate the n hexagons into 
N classes such that hexagons of a particular class are plane 
uncovered areas and the ones of other iV — 1 classes are occu- 
pied by opaque obstacles. We need a class arrangement such 
that no two plane hexagons are adjacent to each other, since a 
single sensor would be sufficient to cover both in that case. The 
problem is the same as the frequency allocation problem to a 
location partitioned into hexagonal cells where we ensure that 
adjacent cells do not share the same frequency. From standard 
results on the problem, we know that N — i'^+ij+j'^, i, j £ 
where i and j are the number of hexagons traversed along two 
axes separated by 120°, to reach a cell of the same frequency 
as the starting cell. We also know that the densest packing of 
cochannel cells occurs for = 3. 

Maximum number of sensors required 

n 

= max — 

n 
~ 3 

The lower bound is attained if each of the anomalous hexagon 
requires no more than a single sensor for its coverage. ■ 

D. Algorithm and Convergence 

The algorithm proposed for the opaque obstacle problem 
has the same ingredients as the one stated for coverage with 
transparent obstacles. However, some extra computation arises 
due to the need to compute the RSP of each point involved in 
a particular placement. As before, the entire area is tessellated 
with hexagons and the anomalous hexagons are identified. 
A 'cluster' is defined as a group of anomalous hexagons 
contiguous to one another. In a particular cluster, let A be the 
area interspersed with opaque obstacles and requiring coverage 
and A be its area. Let M be the number of clusters after the 
initial hexagonal placement and iV be the number of lattice 
points within each hexagon. Let Ui (A) be the area remaining 
to be covered in a region with area A after placement of a 
tessellation whose centres are moved in accordance with the 




shift I E L'. Let us define X to be the set of all coordinates of 
centres of hexagons which result from the application of the 
shift / on the primary tessellation. Let x be an element of X 
and Hr^ {x) refers to the hexagon of circumradius Vs centred 
at X. 



UiiA)^A-A 



[j{H{x%eA)f]RSP{x))f]A 



X 



(8) 



The equation is nearly the same as for transparent obstacles, 
except for the RSP factor coming into play. The factor 
accounts for the fact that the area that gets covered in a 
particular iteration must lie within the RSP of the sensor, 
besides satisfying the usual criteria. Figure Q compares the 
performance of the proposed algorithm with an identical input 
polygon in [12|. The superiority of the proposed algorithm 
is clearly evident. The convergence results derived in the 
previous section also hold for this algorithm, with the sole 
difference lying in the definition of Ui in Equation [8] 

VI. Conclusion and Future Work 

The results presented in this paper also lend themselves 
to solving the problem of fc-coverage, for k=2, 3. Figure [S] 
depicts the standard shift-based procedure for achieving 2 and 
3-coverage. It is to be noted that the procedure ensures that the 
sensors of different tessellations are spaced sufficiently apart. 

2 and 3-coverage can be achieved in the presence of 
obstacles by ensuring that the shifted tessellations are all 1- 
covered, which is in turn accomplished by the application 
of the proposed algorithm separately on each of the shifted 
tessellations. A recent work |15| gives asymptotically optimal 
solutions for A;-coverage problems of arbitrary k for plane 
areas. A possible future extension of the work of this paper 
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